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The Thouless conjecture states that the average conductance of a disordered metallic 
sample in the diffusive regime can be related to the sensitivity of the sample's spectrum to 
a change in the boundary conditions. 

Here we present results of a direct numerical study of the conjecture for the Anderson 
model. They were obtained by calculating the Landauer-Biittiker conductance gi for a sam- 
ple connected to perfect leads and the distribution of level curvatures for the same sample 
in an isolated ring geometry, when the ring is pierced by an Aharonov-Bohm flux. In the 
diffusive regime (L ^> Z e ) the average conductance (gi) is proportional to the mean absolute 
curvature (|c|): (gi) = 7r(|c|)/A, provided the system size L is large enough, so that the 
contact resistance can be neglected. l e is the elastic mean free path, A is the mean level 
spacing. When approaching the ballistic regime, the limitation of the conductance due to 
the contact resistance becomes essential and expresses itself in a deviation from the above 
proportionality. However, in both regimes and for all system sizes the same proportionality 
is recovered when the contact resistance is subtracted from the inverse conductance, show- 
ing that the curvatures measure the conductance in the bulk. In the localized regime, the 
mean logarithm of the absolute curvature and the mean logarithm of the Landauer-Biittiker 
conductance are proportional. 



I. INTRODUCTION 



It has been shown by Thouless and Edwards in the 70 's that the conductance of a dis- 
ordered diffusive system can be related to the dependence of the energy levels to a change 
in the boundary conditions 0,0] ■ A physical realization of this change in the boundary 
conditions is made when the sample is closed to a ring and pierced by a Aharonov-Bohm 
flux (f). In such a case the wave function must obey the condition ip(x + L) = if)(x)e tv where 
1] = 27r0/0o- 0o = h/e is the flux quantum. Thouless found that the average conductance is 
proportional to the width of the distribution of the curvatures of the energy levels when r\ 
is the perturbation parameter. This relation is based on the similar structure of the Kubo 
expression for the conductance and of the curvature of energy levels when the boundary 
conditions are changed. In the diffusive regime, the average absolute curvature of the en- 
ergy levels is proportional to the diffusion coefficient and thus to the conductance. Although 
this conjecture was derived under oversimplified assumptions, its basic idea has proven to 
be very powerful and has become a keystone in our understanding of localization. 

Meanwhile, other measures of the conductance in terms of the response of the system's 
energy spectrum to a change of the boundary conditions have been derived. These spectral 
measures of the conductance were all obtained for the diffusive regime. They all measure 
the diffusion constant in the bulk, and are in this sense equivalent to each other. Indeed, 
it has been shown that the Kubo conductance can be expressed as the average square of 
the first derivative of the energy levels with respect to flux f||-|5[]. The average has to be 
taken over both flux and disorder. The relation between these quantities and the curvature 
distribution was studied both numerically and analytically |4],|6|,|7j] . 

However, so far very little effort has been devoted to comparing spectral measures of the 
conductance to the customary conductance formulas based on transport considerations, like 
in particular the Landauer-Biittiker formula. Exceptions are Ref. ||, where small ID sys- 
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terns were treated numerically, and Ref. || devoted to the study of random banded matrices. 
This is surprising, since spectral measures of the conductance have been widely used in the 
mesoscopic community. Yet, it is well known that the conductance of a mesoscopic sample 
depends sensitively on the measurement geometry and on the way the leads are attached 
to the sample. On the other hand, the spectral measures of conductance are completely 
insensitive in this respect, since the sample is closed to a ring and no leads are attached. 
The system is then in fact a different one. Its spectrum is discrete when the system is finite, 
whereas the system attached to the leads always has a continuous spectrum. 



How can it then be possible that the width of the curvature distribution measures the 
conductance obtained in a transport measurement? In this paper we find that the appro- 
priate conductance is obtained in a situation where the system is connected to leads with 
the same transverse width and with maximal transmission coefficient. We refer to this sit- 
uation as "maximal coupling". It is analogous to the "matching wire" condition defined 
by Economou and Soukoulis |T7| and also used in ref. [[J. In this case the mean absolute 
curvature is proportional to the Landauer-Buttiker conductance. This proportionality holds 
if the system is large enough and the disorder strong enough such that the contact resistance 
can be neglected. In the ballistic regime, the contact resistance is always important and de- 
stroys the proportionality. However, after subtracting the contact resistance from the total 
inverse conductance, the remaining bulk conductance is proportional to the mean absolute 
curvature. The proportionality coefficient is the same in the diffusive and in the ballistic 
regime for all system sizes, as long as the samples have similar geometry. This shows that 
the curvatures measure a bulk conductance. In the localized regime, it is the mean logarithm 
of the absolute curvatures and the conductances that are proportional to each other. 



In the next section we briefly review the different definitions of the conductance. In sec- 
tion JTJ we discuss our numerical method and in section [TV] we present the results of extensive 
numerical simulations, in which we calculated the level curvatures from a perturbative for- 
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inula and the conductances with the Landauer-Biittiker formula. The conductances extend 
over seven orders of magnitude and allow us to study the diffusive as well as the ballistic 
and localized regimes. We conclude in section |V|. 



II. KUBO, LANDAUER AND THOULESS CONDUCTANCES 

A. Kubo and Thouless conductances 

Let us first recall the main line of the Thouless derivation. On the one hand, the d.c. 
Kubo conductivity can be written as [Kj : 



° = -5x7 E \P«p\ 2 6(E F - e a )6(E F - ep ) (1) 

where e a are single energy levels and p a p = (a\p\/3) are the matrix element of the momentum 
operator. V = L d is the volume and Ep the Fermi energy. Strictly speaking, this expression 



is zero for a finite system | H| . To get a finite a, the 8 functions must have a finite width 
larger than the inter-level spacing. Under this condition and assuming that the matrix 
elements p a p are decorrelated from the e a JT2], the average conductivity is given by: 



ne 2 hL 2d ~ 2 
m 2 



{G K ) = (a)L d ~ 2 = " (\p af} \ 2 )pl . (2) 



(...) represents an average over the disorder. p is the average density of states per unit 
volume at the Fermi energy. The dimensionless conductance (gx) can be written as fT3| : 

(Gk) E, 



k) = ^t = 2^ (3) 



e 2 /h A 

where A = 1/ (poL d ) is the mean level spacing, and the Thouless energy E c is given by 
E c = hD/L 2 , D being the diffusion coefficient [0. The second equality in formula @ is 
nothing but the Einstein relation a = e 2 Dp Q . All these quantities are defined for a given 
Fermi energy and will in general depend on Ep. 

On the other hand, under the change in the boundary conditions tp(x + L) = if)(x)e lv , the 
curvature of a given energy level e a at the origin (77 = 0) is given exactly by perturbation 
expansion in rj: 



fd 2 e a \ = _^ + \p^l 
° a V dr l 2 ) v =o mL2 m2L2 ^a £ *- £ p' 

In order to relate the width of the curvature distribution to the diffusion coefficient, Thouless 

assumed first that the energy levels e a are not correlated with the matrix elements p a p. 

Replacing then \p a /3\ 2 by its average value, the distribution of the curvatures is that of 

l/(e a — £/?)■ Secondly, assuming that the energy levels themselves are not correlated, the 

sum in eq.(^) gives rise to a Levy law for the distribution of the curvatures in the limit of 

infinitely many levels |T|. It has the Cauchy form P(c) = (7o/ 7r )/(7o + ° 2 ) with a width 70 

given by 

^° = ~ 2T2 — x — • ( b > 

Comparison between the equations (§) and @ gives the relation between the dimensionless 
average conductance (gx) = (GK)h/e 2 and the width of the distribution of curvatures, 
known as the Thouless relation @: (gx) = 7r^. However, it is now known that the energy 
levels are strongly correlated in a metal so that the curvature distribution does not have the 
Cauchy form. Instead, it is given by: 

P P( C ) = ( 7 | + c 2)(/3+2)/2 ( 6 ) 

Here, (5 = 1 if there is time-reversal symmetry and (3 = 2 if time-reversal symmetry is 
broken. Np is a normalization coefficient. This form has first been guessed by Zakrzewski 
and Delande to fit numerical calculations on various models exhibiting chaotic spec- 



tra. It has been proven analytically by von Oppen |15[] for random matrices of the form 



H(X) = H + XK where H and K are random matrices belonging to the same symmetry 
class {(3 = 1 for the Gaussian Orthogonal Ensemble, GOE; (3 = 2 for the Gaussian Uni- 
tary Ensemble, GUE). A is the perturbation parameter. Recent numerical calculations have 
shown that this distribution is also characteristic of metallic spectra when the perturbation 
parameter is an AB flux @. In particular, in the limit where — > 0, the distribution 
is still the GOE distribution (/3 = 1 in eq. (ffl)) ||. This has been proven analytically by 



Fyodorov and Sommers who also found that there are no corrections of order A/E c 



The normalized distribution in zero field is thus: 

Fyodorov and Sommers have shown that the width of this distribution can be related to the 
diffusion coefficient [J/J. They find that the width 71 of this distribution for a three dimen- 
sional ring is given by 71 = 2E C . Using now eq. (Q) relating E c to the Kubo conductance, 
one deduces 71 = A(g K )/n. To characterize this width it is convenient to introduce the 
average absolute curvature (|c|) = 71, so that the Thouless conductance defined as 

(9t) = ttM (8) 

equals the Kubo conductance: (gx) = (9k)- 

It has to be noted: first, all conductances considered so far are average conductances, the 
average extending over disorder realizations. Secondly, the equation (g?) = (gx) holds so 
far only for the diffusive regime, since 71 = 2E C was derived in |?J for the diffusive regime. 
We will see that the equation has to be modified in the ballistic and localized regimes. 

B. Landauer-Biittiker conductances 

Another way to express the conductance has been introduced by Landauer. He related 
this quantity to the scattering properties of the disordered system, when it is connected 
to incoherent reservoirs through ideal leads. This approach ideally suits transport through 
finite mesoscopic systems and shows the importance of the measurement geometry. For one 



dimension, Landauer derived the dimensionless conductance gi [16] 



9L = ^Jh = T-T- (9) 

This conductance is the ratio Gl = I / '(/M — where ha and fiB are the chemical potentials 
of ideal leads attached to the barrier. T is the transmission coefficient through the disordered 
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region, diverges for an ideal, clean sample. On the other hand, Economou and Soukoulis, 



trying to derive this formula from linear response theory (Kubo formula), found |I7|: 



instead. In this case, Gl = I/i^i ~ ^2)1 where fi\ and \xi are the chemical potentials of the 
reservoirs Eq.@ describes a four-terminal measurement in one dimension, that is a 



measurement with separate current and voltage probes ||18|| . gt = T describes a two-probe 



measurement, where only two leads are attached to the sample and serve as current and 
voltage probes at the same time. The remaining finite resistance at zero disorder (gi = 1) 
is a "contact resistance" which has its origin in the coupling of the sample to the incoherent 



reservoirs 



[p!T| , |19[ . This resistance cannot be avoided in a two-probe measurement. One may 
therefore think of the total resistance Gj} in such a ID two-probe geometry as being the 
sum of the contact resistance h/e 2 and a "bulk resistance" Gj} . The latter vanishes when 
the disorder goes to zero and is identical with the original Landauer contribution: 



h 



Gl l = G L l + - (11) 



Fisher and Lee generalized eq.(|Il]) to the multi-channel case plf : 



M 



g L = 9K = T, T i = tTtt+ - ( 12 ) 



i=l 



Ti is the total transmission probability in the i th channel, t the transmission matrix and M 
the number of channels. When the disorder in the sample goes to zero, gi is limited by 



the number of open channels. Today there is a general agreement that eq.(|12|) describes a 
two-probe measurement in a multi-channel geometry. In our numerical simulations we will 
focus on this situation and use eq.(|l2D for the numerical evaluation of the conductance. 
Comparing equations (||) and (p!2|), one gets a relation between the average Landauer- 
Biittiker conductance and the width of the curvature distribution: 

<^> = 7rM (13) 
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In the two-probe multi-channel geometry we will consider in the following one might again 
decompose the total resistance into a sum of a contact resistance plus a bulk resistance, the 
latter being entirely due to the motion in the bulk of the sample. In straight generalization 
of eq. (J!]]) , it is then natural to define the bulk conductance Gl = j[9l by 

Gl x = Gl 1 +R c , (14) 

where R c = h/(Me 2 ) is the "contact resistance" for the multi-channel system [11,21 1. We 
have then 

~ 9L = 7^YK- (15) 

1 M 

In the diffusive regime, the effective number of conducting channels, M e ff = J2Ti, is much 



smaller than M: M e ff = Ml e /L where l e is the elastic mean free path p3| |. Consequently, 



gi and gi are almost identical in the diffusive regime, the relative deviations being of order 
l e /L. However, in the ballistic regime they behave very differently: cjl — > oo and — > M 
in the limit of zero disorder. 

For more than one channel gi has not the simple and general interpretation of the 
conductance measured in a four-probe measurement. Indeed, that in the multi-channel case 
not only the number of leads but also the way (e.g. under what angles) they are attached 
influences the measured conductance, such that a general four-probe formula might not even 
exist |p2| . Similarly (|c|) cannot correspond to any particular four-probe conductances, since 
it is an intrinsic property of the disordered region. We will show that for a finite system 
(|c|) is proportional to {gi)'- 

i~ \ d c D tit* 

m) (16) 

III. THE NUMERICAL METHOD 



The starting point of our analysis is the Anderson tight-binding Hamiltonian if of a 
disordered system on a square lattice of L x x L y x L z sites. For the curvature calculation 
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the system is closed to a ring and pierced by an Aharonov-Bohm flux <p: 

# = £e,|z>(*|+u$: \i)(j\+u £ (e iv \i)(j\ + h.c). (17) 

i <ij> <ij> 

i x =L,j x =l 

The Ci are distributed uniformly and independently in an interval between —w/2 and w/2. 
< ij > denote next nearest neighbors, u is the hopping matrix element which we set equal 
to one in the following, and w is the disorder parameter. The last sum in eq.(|l7]) is over 
the set of sites on the two boundaries limiting the open sample in x-direction. Hopping 
between these boundary sites arises when the system is closed to a ring and includes a phase 
factor e ir) . For = or entire multiples of the flux quantum, one recovers periodic boundary 
conditions. 

For the calculation of gL, the system is open and coupled to perfect leads. The last sum 
in eq.(|lj]) is then missing. This is the only difference between the two Hamiltonians. In 
particular, for the numerical implementation the same random number generator was used 
for the diagonal matrix elements in both situations. 



A. Curvatures 

In the diffusive regime, the curvatures can be evaluated by replacing differentials by small 
flux differences whose values are varied for control in a suitable way || . This procedure has 
the numerical advantage that only eigenvalues, not the eigenvectors are needed. However, it 
is very difficult to control in the ballistic and in the localized regime. We adopted therefore 
a routine based on an exact perturbative formula corresponding to eq.([|). In fact, treating 
rj in eq. ( |17D as a perturbation up to second order, one finds for the curvatures at zero flux 

f V 



c„ 1 

12 ( e a\i)(j\e a ) + Y 12 (ep\i)(j\e a ) - (e^\j)(i\e a ) 



;is) 



l ix = L x jjx—^- \ix = Lx ijx — 1 / 

where e a and \e a ) denote the eigenvalues and eigenvectors of the Hamiltonian at zero flux, 
respectively. Higher order terms vanish since 77 = 0. In the two directions perpendicular 
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to the transport direction, periodic boundary conditions were used. Formula fll8|) is exact 
as long as ^ e a . Thus, for a finite system, where level repulsion is always present at 
sufficiently small energy scales |[25|| , (|T8| ) remains valid also in the localized and the ballistic 
regime. Besides rounding errors which can be neglected here the only remaining errors in 
the calculation of (|c Q |) are statistical errors that can be controlled by increasing the number 
of disorder realizations. We used up to 1000 disorder realizations for system sizes of 6 x 6 x 6 
sites and still about hundred for 10 x 10 x 10 sites. Relatively, the remaining statistical errors 
in the diffusive and ballistic regimes were of the order of 10 -2 , which we checked by varying 
the number of disorder realizations. As eq . fll8|) indicates, all eigenvalues and eigenvectors 
are needed for the calculation of a single curvature. Realizations where our Lanczos routine 
failed to find all eigenvalues and eigenvectors were therefore discarded. 



B. Conductances 

The conductance Ql was calculated from eq. (|T2|) by the Green's function recursion tech- 
nique The Green's function connecting the 2 ends of a strip can be calculated recursively 
using the equations 

(19) 

r (N) _ r (N-l) r (N) , 9m 

where G^ N represents the sub-matrix of the Green's function between sites on the iVth slice 
of a strip of length N, G[ N ^ is the corresponding sub-matrix between sites on the 1st and 
iVth slices, and represents the Hamiltonian of the iVth slice alone. The system can be 
embedded in semi-infinite leads by choosing the initial values of the 2 Green's functions to 
represent the end of a semi-infinite wire and by adding a final slice for which the Hamiltonian 
of the slice is replaced by the self-energy matrix for another semi-infinite wire. Having the 
Green's functions one can derive the transmission matrix t and then the conductance 
9l- 



G 



(AT) 
N.N 



Z - Hjv - u+G^/l^u 
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IV. RESULTS 



A. Energy dependence 

Without averaging over energy, both (g^) and (gi) are energy dependent: (gx^E)) and 
(gz,(E)). The variation of (g^E)) is smooth and is due to the energy dependence of the DOS 
and of the number of channels M. The energy dependence of {gT{E)) arises from the varia- 
tion with energy of both (|c|) and A(E), where the latter quantity is the mean level spacing 
at a given energy (averaged over disorder only). In order to get the conductance at a given 
energy, we therefore rescaled the curvatures with an energy dependent A: c a = c Q /A(e a ). 
The disorder averaged DOS, 1/A(E) , was obtained by the standard method of fitting the 
spectral staircase (integrated DOS) to a polynomial. 

After averaging over 1000 disorder realizations (in the case of systems with 6x6x6 
sites), the fluctuations of (|c a |) = (|c a |)/A( function of energy turned out to be 

still much more pronounced than those of (gi{E)). This is not too surprising, as it is 
well known that in the diffusive regime the conductance distribution (which is a universal 
Gaussian distribution [^,^] with a width of the order of the conductance quantum) and 
the curvature distribution (see section |) are very different. Thus, when using just one 
disorder realization, the fluctuations of the function |c a (e a )| will be much larger than those 
of (g(E)), due to the long 1/c^ tails of the curvature distribution. For a finite number of 
realizations this difference will still persist, and only when averaging over infinitely many 
disorder realizations the energy dependence of (gi) should follow that of (|c a |). Having 
in mind that even 1000 disorder realizations did not suffice to reduce the fluctuations of 
{gT{E)) to a level comparable to those of {gi(E)), it seems very difficult to check the 
Thouless conjecture in the stronger sense for a given energy with the current computing 
power available. We therefore averaged gi and gx not only over the realizations but over a 
band of energy AE comprising typically about the central half of the spectrum as well. We 
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checked that increasing the size of the system or the number of realizations allow to decrease 
AE to obtain the same results. This suggests that our results are independent of AE. In 
the following, (. . .) will stand for the combined disorder and energy average. Care was taken 
in order to average both curvatures and conductances over exactly the same energy interval. 



B. Curvature distribution 

In the diffusive regime, the distribution of the curvatures is well described by eq. ([?]). 
Thus, (|c|) is a good measure of the width of the curvature distribution. 
Outside the diffusive regime the curvature distribution was not known so far, and one might 
wonder whether (|c|) is still well defined. We therefore calculated P(c) numerically for both 
the ballistic and the localized regimes. Fig.[l] shows P(c) for a system in the ballistic regime 
(6x6x6 sites, w = 1.0, 4000 disorder realizations) and the prediction of eq.(^), where 71 
was determined as 71 = (|c|) (no fitting parameter). Eq.(|7]) works well for large curvatures 
and shows that in the ballistic regime P(c) has 1/c 3 tails as in the diffusive regime. For 
small curvatures deviations from eq.(|7]) in the form of non-universal features appear and the 
distribution develops two maxima. A relative minimum appears at zero curvature. These 
deviations become even more pronounced for smaller disorder. Altogether we conclude that 
(|c|) can still serve as a measure for the width of the curvature distribution, even in the 
ballistic regime. 



In the localized regime at least two different numerical works favor a log-normal curvature 



distribution [pq , p^| . Analytical evidence for a log-normal distribution at least for small 
curvatures in ID is given in ||30|| . On the other hand, one might suspect that Thouless' 
original result of a Cauchy distribution due to uncorrelated eigenvalues might apply to the 
localized regime. Such a distribution would of course spoil the use of (|c|) as a measure of 
the Thouless conductance. We therefore reexamined this question numerically. As shown 
in Fig.0, a Cauchy distribution can be ruled out: for large curvatures, the distribution falls 



12 



off faster than 1/c 2 . This can probably be explained by the fact that the eigenvalues and 
the eigenvectors are strongly correlated for large disorder in contrast to what was assumed 
by Thouless in the derivation of his formula. On the other hand, Fig.^ shows that a log- 
normal distribution does not fit perfectly either. Rather large deviations are visible for large 
curvatures. We will address this question in more detail in a future work. Nevertheless, 
we can conclude from Fig.^| that both (|c|) and (In |c|) are well-defined quantities in the 
localized regime [311. 



C. Disorder Dependence 

Before discussing the disorder regimes separately, we display in Fig.^ an overall plot of 
the disorder dependence of (\c a \) and («7l). Several points can be observed immediately: 
First of all, (|c a |) diverges for small disorder like 1/w 2 . This is a well-known fact which can 



be derived from perturbation theory (first Born approximation [3~2||). Also, (\c a \) has the 
right scaling behavior of a conductance. In 3D, in the ballistic and diffusive regimes, (\c a \) 
increases proportionally to the system size L within the parameter range provided (L = 6 
to L = 10). In the localized regime it decays with the system size. In 3D there is a critical 
value w c ~ 16.5 where (\c a \) becomes independent of the system size, thus indicating the 
position of the metal-insulator transition (MIT). Within the error bars it coincides with the 
well-known value found by MacKinnon and Kramer, who examined the scaling behavior 



of the transmission through disordered samples |33]. We also checked that in 2D (\c a \) is 



independent of L in the diffusive and ballistic regimes. 



1. Ballistic and Diffusive Regimes 

Fig.| shows that (g^) obeys the same scaling behavior as (|c a |). However, the disorder 
dependence of (gi) and (\c a \) is rather different. Even in the diffusive regime, (gi) follows 
(|c a |) only over a small disorder interval close to the metal-insulator transition. The inter- 
val's width increases with the system size, but for all system sizes the discrepancy becomes 
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very pronounced in the ballistic regime, where (gi) converges to a constant value, whereas 
(\c a \) keeps diverging. 



Following our discussion of section [H] this result is not surprising. In the ballistic regime, 
D formally diverges, as does ( | Cq, |) . Any limitation of the conductance due to the coupling of 
the sample to the environment must then result in a deviation from the conjectured propor- 
tionality between (\c a \) and (gi)- Clearly, the discrepancy in the lower disorder limit of the 
diffusive regime is already caused by the cross-over to constant {gi) due to the boundary 
resistance. 



In order to improve the agreement of (\c a \) with the conductance, the latter has to be 
defined such that it does not incorporate the contact resistance. We therefore also compared 
the disorder dependence of (gi,) with the one of (|c a |). As explained in sec.0, (gi) does not 
contain the contribution of the boundary resistance and should be a measure of the bulk 
conductance. It will therefore also diverge when the disorder vanishes. Whereas it is not 
clear from the beginning that this divergence will be of the same kind as the one of (|c Q |), 
Fig. |5] shows that (gi) diverges for small w indeed with the same power as (|c Q |). Both 
curves follow each other from the diffusive regime until far into the ballistic regime. 



In Fig. | we have plotted (gi) as function of (|cq,|). The MIT is given in this plot by the 
point where 47r 2 (|c Q |) ~ 4.1. The points from all sample-sizes considered in 3D now fall on 
one straight line with slope one and this in the diffusive as well as in the ballistic regime 



34j . A fit to a linear law gives 

(g L ) = (0.99 ± 0.04)7r(|c a |) - 0.029 ± 0.008 (21) 

in remarkable agreement with eq.(^) |[36| . The error bars were obtained as standard devia- 
tions from the three system sizes considered. The conclusion is therefore that (\c a \) measures 
the bulk conductance (gi) in both the diffusive and the ballistic regime. 
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In a recent paper, Casati et al. || also study the relation between Landauer conductance 
and curvature distributions for band random matrices and they find the relation: 

(g L ) = (7.5 ± A)K W (22) 

where they define K av as the geometric average exp(ln(|c|/A)) ||. For a distribution of 
curvatures like [7], the geometric average is related to the arithmetic average by 



K av = exp(ln^) = i(^) 



Using the relation [13] we get the following result: 

(g L ) = 2nJC av 

This factor 2tt satisfactorily explains the numerical result {eq. ^} found by these authors. 

2. Localized regime 

Fig.[] shows that the power between (\c a \) and (gi) changes at the MIT. We obtain 
approximately (gi) oc (|c a |) L2 . However, in the localized regime (|c a |) and (gi) might 
not be the right quantities to look at. At least from the conductance it is known that in 
this regime the function with the right scaling behavior is (higi) ||37|| , not {ql}- Since the 



favored log-normal distribution of curvatures is due to the same reason as the log-normal 
distribution of the conductances, namely the exponentially decaying wave functions with 
normally distributed localization length, one might suspect that a similar statement holds 
for the curvatures as well. We therefore also examined (lng^) and (In |c| ) as functions of 
disorder and system size. Fig.[7| shows the result: first of all both quantities are proportional 
to the system size, and secondly proportional to each other. Plotting (lng^) versus (In | c|) 
(see Fig.|) yields a straight line that is best approximated by the linear law 

(ln# L ) ~1.7(ln|c|)-2.5. (23) 

Again, the validity of this equation extends over several orders of magnitude of the conduc- 
tance. However, the prefactor seems to decrease slightly but systematically with the system 
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size. A law like ^3] was also reported in || with a similar prefactor (1.73) for banded random 
matrices. Nevertheless, Ref. || also reports a prefactor 2.0 for an Anderson model. We do 
not have any explanation for this difference besides the fact that in contrast to our box- 
distributed disorder the disorder was Gaussian distributed in ||. The found behavior is a 
priori surprising. Assuming that in the localized regime the flux dependence of each energy 



level is purely sinusoidal |38| one deduces that (i 2 ((p)) oc (c 2 ), where i a (<f) = —de a /dip and 



the overline indicates a flux average. Since (i 2 ((p)} can be related to the Kubo conductance, 
one would expect a quadratic relation (gi) oc (c 2 ) |39 . 



V. CONCLUSION 

We have examined numerically the relation between level curvatures and conductances 
for disordered systems. We showed that in the diffusive regime a proportionality between 
the dimensionless mean absolute curvature (\c a \) and the average Landauer-Biittiker con- 
ductance (gi) holds if the system is large enough so that the influence of the boundary 
resistance can be neglected. In the ballistic regime, the boundary resistance can never be 
neglected and leads to a strong violation of the proportionality. In the limit of zero disorder 
it completely dominates the total resistance and limits (gi) to the number of open channels, 
whereas (|c a |) diverges in the same limit. However, for all system sizes a proportionality 
between a properly defined bulk conductance {gi) and (\c a \) could be established that holds 
in the same form in the diffusive and ballistic regimes. This shows that in these regimes level 
curvatures measure a conductance that is entirely due to the dynamics in the bulk of the 
sample and therefore not influenced by details of the measurement setup, like the number 
of leads and the way they are attached. In the localized regime, we found a proportionality 
between (In | c|) and (lng^). 
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FIGURES 

FIG. 1. The curvature distribution in the ballistic regime. The full line is the prediction of 
eq.(^) known to be valid in the diffusive regime. Deviations at small curvature are visible (see 
inset). 6x6x6 sites, w = 1.0, 4000 disorder realizations. 

FIG. 2. The curvature distribution in the localized regime with a fit to a Cauchy distribution 
(dashed line). This plot shows that P(c) decays faster than 1/c 2 for large curvatures. (6x6x6, 
w = 50, 4000 disorder realizations). 

FIG. 3. Distribution of In |c| in the localized regime with a fit to a Gaussian distribution 
(corresponding to a log-normal distribution for \c\). Same parameters as in Figj2|. 

FIG. 4. The overall disorder dependence of (\c a \) (top) and {g£) (bottom) for different system 
sizes: diamonds 6x6x6, circles 8x8x8, and triangles 10 x 10 x 10. The straight lines of {\c a (w)\) 
in the logarithmic plot correspond to a l/w 2 divergence for small w. Full lines are guides to the 
eye only. For clarity the (\c a {w)\) curves were shifted by an arbitrary factor 4tt 2 . 

FIG. 5. The overall disorder dependence of (\c a \) (top) and (g£) (bottom) for different system 
sizes (same symbols as in FIG.f|). {c)l{w)) diverges in the same manner as (\c a (w)\) for small w. 
Full lines are guides to the eye only. 

FIG. 6. The conductance {g) plotted against (\c a \) for different system sizes (same symbols 
as m Fig. £§. The diffusive regime starts with the critical mean curvature 47r 2 (|c Q |) ~ 4.1. In 
this regime and the ballistic regime the dependence is very well fitted by the same linear law 
(g L ) = (0.99 ± 0.04)vr(|5 Q I) - 0.029 ± 0.008 (full line). 

FIG. 7. (ln|c|) (top) and (Ingi) (bottom) as a function of disorder for three different system 
sizes in the localized regime (same symbols as in Fig. ||). 

FIG. 8. (lng^) as a function of (In |c|) in the localized regime for three different system sizes 
(same symbols as in Fig. |j). 
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